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Abstract 

We discuss calculations of the inclusive total cross section for heavy quark pro¬ 
duction at hadron collider energies within the context of perturbative quan¬ 
tum chromodynamics, including resummation of the effects of initial-state soft 
gluon radiation to all orders in the strong coupling strength. We resum the 
universal leading-logarithm contributions, and we restrict our integrations to 
the region of phase space that is demonstrably perturbative. We include a 
detailed comparison of the differences between ours and other methods. We 
provide predictions of the physical cross section as a function of the heavy 
quark mass in proton-antiproton reactions at center-of-mass energies of 1.8 
and 2.0 TeV, and we discuss estimated uncertainties. 
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I. INTRODUCTION AND MOTIVATION 


In this report we present and discuss calculations carried out in perturbative quantum 
chromodynamics (QCD) of the inclusive cross section for the production of heavy quark- 
antiquark pairs in hadron reactions JIJ-Q. In most of the paper, we identify the heavy 
quark as a top quark, t, but the results are valid more generally as long at the mass of 
the quark is sufficiently heavy. For example, they should apply as well to production of a 
fourth-generation quark, such as a postulated b'. 

In inclusive hadron interactions at collider energies, hi + h 2 —> t+t+X, ti pair production 
proceeds through partonic hard-scattering processes involving initial-state light quarks q and 
gluons g. In lowest-order QCD, at the two partonic subprocesses are q + q —>• t + i and 

g + g ^ t + t. Calculations of the cross section through next-to-leading order, 0{a^), involve 
gluonic radiative corrections to these lowest-order subprocesses as well as contributions from 
the q + g initial state 0. A complete hxed-order calculation at order C>(a"),n > 4 does not 
exist. In this paper, we do not examine mechanisms for the production of single top quarks 
or antiquarks p. 

The physical cross section for each production channel is obtained through the factor¬ 
ization theorem. 




r-1 


drj^ 




Aw? 


(1+h),h" 




( 1 ) 


The square of the total hadronic center-of-mass energy is S', and the square of the partonic 
center-of-mass energy is s. The mass of the heavy quark is m, and /i is the common 


renormalization/factorization scale of the problem. The variable g = 


4m^ 


— 1 measures the 


distance from the partonic threshold. The indices ij G {qq, gg} denote the initial parton 
channel. The partonic cross section aij{g,m‘^, /j?) is obtained commonly from fixed-order 
QCD calculations P, or, as described here, from calculations that go beyond fixed-order 
perturbation theory through the inclusion of gluon resummation pl-ffl. The parton flux 


is = fy where n'^) is the density of partons of 
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type i in hadron hi. We use the notation a(/i) = as(/i)/7r. Unless otherwise specihed, 
a = = m) throughout this paper. The total physical cross section is obtained after 

incoherent addition of the contributions from the the qq and gg production channels. In this 
paper, we ignore the small contribution from the qg channel. 

Comparison of the partonic cross section at next-to-leading order with its lowest-order 
value reveals that the ratio becomes very large in the near-threshold region. Indeed, as 
r; —0, the “i^-factor” at the partonic level K{rj) grows in proportion to a\T?{7]). An 
illustration of this behavior may be seen in Fig. 7 of Ref. [3]. The very large mass of the top 
quark notwithstanding, the large ratio K[rj] makes it evident that the next-to-leading order 
result does not necessarily provide a trustworthy quantitative prediction of the top quark 
production cross section at the energy of the Tevatron collider. The large ratio casts doubt 
on the reliability of simple fixed-order perturbation theory for physical processes for which 
the near-threshold region in the subenergy variable contributes significantly to the physical 
cross section. Top quark production at the Fermilab Tevatron is one such process, because 
the top mass is relatively large compared to the energy available. Other examples include 
the production of hadronic jets that carry large values of transverse momentum and the 
production of pairs of supersymmetric particles with large mass. To obtain more dependable 
theoretical estimates of the cross section in perturbative QCD, it is important hrst to identify 
and isolate the terms that provide the large next-to-leading order enhancement and then to 
resum these effects to all orders in the strong coupling strength. 

We begin in Sec. II with the motivation for the inclusion of the effects of intial state 
soft gluon radiation to all orders in the QCD coupling strength, and we review the general 
formalism of resummation. In Sec. Ill, we outline the method and domain of applicability of 
perturbative resummation that we developed in the past year [0,Q. We present predictions 
in Sec. IV of the physical cross section as a function of the heavy quark mass in proton- 
antiproton reactions at center-of-mass energies of 1.8 and 2.0 TeV, and we discuss estimated 
uncertainties. Our calculation is in good agreement with the measured cross section at the 
reported mass of the top quark [^. At m = 175 GeV and = 1.8 TeV, the all-orders 
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resummed cross section is about 9% greater than the next-to-leading order value. Since the 
large threshold logarithms are mastered by resummation, the theoretical reliability of the 
resummed result is considerably greater than that of the hxed order calculation. At other 
values of m and y/S, where the ratio m/'/S is larger, the numerical effects of resummation 
can be more signihcant. In Secs. V and VI, we compare our approach and results with other 
methods [|l|,0 and address criticisms that have been made [§]. The difference between our 
approach and that of Ref. 0 resides in the treatment of subleading logarithmic contributions, 
and we explain our reasons for preferring our method. Conclusions are summarized in Sec. 
VII. 


II. GLUON RADIATION AND RESUMMATION 

The origin of the large threshold enhancement may be traced to initial-state gluonic 
radiative corrections to the lowest-order channels. To avoid misunderstanding, we remark 
that we are calculating the inclusive total cross section for the production of a top quark- 
antiquark pair, i.e., the total cross section for t + t + anything. The partonic subenergy 
threshold in question is the threshold for t + t+ any number of gluons. This coincides with 
the threshold in the invariant mass of the t + t system for the lowest order subprocesses only. 

To specify the kinematic variables, we consider the two-to-three parton subprocess i{ki) + 
j{k 2 ) i{pi) + i{p 2 ) + 9 {k). We dehne the variable z through the partonic invariants |l[] 

s= {ki + k 2 Y, ti = {k 2 - P 2 Y - “1 = (^1 “ P'if - (1 - z)m‘^ = s + ti + ui. (2) 

Alternatively, {1 — z) = In the limit that z —>■ 1, the radiated gluon g{k) carries zero 

momentum. After cancellation of soft singularities and factorization of collinear singularities 
in 0{a^), there are left-over integrable logarithmic contributions to the cross section asso¬ 
ciated with initial-state gluon radiation. The contributions of interest here, often expressed 
in terms of “plus” distributions, are proportional to ln(l — z). These logarithmic terms are 
vestiges of the canceled infrared singularities. 
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The partonic cross section may be expressed generally as 



dzHij{z, a)alj{ri, m^, z). 


( 3 ) 


We work in the MS factorization scheme in which the q, q and g densities and the 
next-to-leading order partonic cross sections are dehned unambiguously. The lower limit 
of integration, Zmin = 1 — 4(1 + rj) + 4^1 -|- rj, is set by kinematics. The derivative 
z) = d{alj\ri,m‘^,z))/dz, and d® is the lowest-order 0(a‘^) partonic cross sec¬ 
tion expressed in terms of inelastic kinematic variables to account for the emitted radiation. 

Keeping only the leading logarithmic contributions through 0{a^), we may approximate 
the total partonic cross section as 


rri^) = j dz^l + a2Cij ln^(l - 

^ ZrnAn 

dzn[^j~^^\z,a)aUT],z,m^) , 


Z'mzn 
nl 


( 4 ) 


where Cqq = Cp = 4/3 and Cgg = Ca = 3. As is illustrated in Fig. 1 , the leading logarithmic 
contribution, integrated over the near-threshold region 1 > ^ > 0 , provides an excellent 
approximation to the exact full next-to-leading order physical cross section as a function of 
the heavy quark mass. 

Although a hxed-order 0{a‘^) calculation of ti pair production does not exist, we may 
invoke universality with massive lepton-pair production (//), the Drell-Yan process, to gen¬ 
eralize Eq. (^) to higher order. In the near-threshold region, the hard kernel becomes 


a) ~ 1 -t- 2aCij ln^(l — z) + 


2 Cj ^*^(1 -z) - -Cijb 2 ln^(l - z) 


( 5 ) 


The coefficient 62 = (HC'a ~ 2 ? 7 ,j)/ 12 , and the number of flavors Uf = 5. We note that 
the leading logarithmic contributions in each order of perturbation theory are all posi¬ 
tive in overall sign [Q so that the leading logarithm threshold enhancement keeps build¬ 
ing in magnitude at each hxed order of perturbation theory. The further enhancement of 
the physical cross section produced by the 0{a'^) leading logarithmic terms in the near¬ 
threshold region is shown in Fig. 1 . At m = 175 GeV, we compute the following ratios of 
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the physical cross sections in the leading logarithmic approximation: =1.22, and 


= 1 - 14 . 


The goal of glnon resnmmation is to snm the series in a” ln^”(l — z) to all orders in a in 
order to obtain a more trustworthy prediction. This procedure has been studied extensively 
for the Drell-Yan process |p, and good agreement with data is achieved. In essentially 
all resnmmation procedures, the large logarithmic contributions are exponentiated into a 
function of the QCD running coupling strength, itself evaluated at a variable momentum 
scale that is a measure of the radiated gluon momentum. For example, in the approach of 
Laenen, Smith, and van Neerven (LSvN) [|l|, the resummed partonic cross section is written 


as 




( 6 ) 


where the exponent 


Eij{z^m?) cx Cija{{l — zY^^m^) ln^(l — z). 


(7) 


We note that in Eq. (|^), the strong coupling strength is evaluated at the variable momentum 
scale (1 — 

Different methods of resnmmation differ in theoretically and phenomenologically im¬ 
portant respects. The set of purely leading monomials a"^ln^”(l — z) in djj exponentiates 
directly, with a evaluated at a hxed large scale /i = m, as may be appreciated from a glance 
at Eq. d^). This simple result does not mean that a theory of resnmmation is redundant, 
even if only leading logarithms are to be resummed. Indeed, straightforward replacement 
of the term within the brackets of Eq. (H) with the exponential of a2Cijln'^{l — z) would 
lead to an exponentially divergent integral (and therefore cross section) since the coefficient 
of the logarithm is positive. The naive approach, therefore, fails from the start, and more 
sophisticated resnmmation approaches must be employed, involving scaling and Lorentz- 
transformation properties of the classes of terms to be summed. The more sophisticated 
approaches are not free from problems, however. Formally, if not explicitly in some ap¬ 
proaches, an integral over the radiated gluon momentum z must be done over regions in 
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which z —i> 0. Therefore, one signihcant distinction among methods has to do with how the 
inevitable “non-perturbative” region is handled in each case. Examination of Eqs. and 
(^) shows that an infrared singularity is encountered in the soft-gluon limit z ^ 1: owing to 
the logarithmic behavior of a{q^), a{q‘^) oc ln“^(g^/AQ(^^), a((l — oo as z ^ 1. 

The infrared singularity is a manifestation of non-perturbative physics. In the approach of 
LSvN, this divergence of the integrand at the upper limit of integration necessitates intro¬ 
duction of the undetermined infrared cutoff (IRC) po in Eq. (|), with Aqcd < Ao < The 
cutoff prevents the integration over z from reaching the Landau pole of the QCD running 
coupling constant. The presence of an extra scale spoils the renormalization group proper¬ 
ties of the overall expression. The unfortunate dependence of the resummed cross section 
on this undetermined cutoff is important numerically since it appears in an exponent [|l|. 
Theoretical uncertainties are not easy to evaluate quantitatively in a method that relies on 
an undetermined infrared cutoff. 


III. PERTURBATIVE RESUMMATION 


The method of resummation we employ is based on a perturbative truncation of 
principal-value resummation (PVR). The principal-value method has an important tech¬ 
nical advantage in that it does not require arbitrary infrared cutoffs, as all Landau-pole 
singularities are by-passed by a Cauchy principal-value prescription. Because extra unde¬ 
termined scales are absent, the method also permits an evaluation its perturbative regime 
of applicability, i.e., the region of the gluon radiation phase space where resummed pertur¬ 
bation theory should be valid. 

To illustrate how infrared cutoffs are avoided in the PVR method, it is useful to begin 
with an expression in moment (n) space for the exponent that resums the ln(l — z) terms 
I0| . Factorization and evolution lead directly to exponentiation in moment space: 

1 „ 1 . 1 


E{n, irY) = 


dx 


X 


n—1 


dX 


1 — X J X 

(1-X)2 


■9 


a 


{Xrn?^ 


( 8 ) 
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The function g{a) is calculable perturbatively, but the behavior of Q;(Am^) leads to a di¬ 
vergence of the integrand when > ^qcd- tame the divergence, a cutoff can be 

introduced in the integral over x or directly in momentum space, in the fashion of LSvN. In 
the principal-value redehnition of resummation, the singularity is avoided by replacement 
of the integral over the real axis x in Eq. by an integral in the complex plane along a 
contour P that has the same endpoints and is symmetric under reflections across the real 


axis: 


{n,m^) = — I d( 


X 


n—1 


dX 


l-C ^ A 


■9 


a 




(9) 


The function is hnite since the Landau pole singularity is by-passed. Moreover, 

\im.n^oo {n,'ni?) = —oo, and, therefore, the corresponding partonic cross section is hnite 
as 2 ; —>• 1 (n —>• -|-cx)). In Eq. (^), all large soft-gluon threshold contributions are included 
through the two-loop running of a. 

Equations (|) and have identical perturbative content, but they have different non- 
perturbative content since the infrared region is treated differently in the two cases. The 
non-perturbative content is not a prediction of perturbative QCD. In our study of top quark 
production, we choose to use the exponent only in the region of phase space in which the 
perturbative content dominates. 

We use the attractive hniteness of Eq. (|^) to derive a perturbative asymptotic represen¬ 
tation of E{x, a{m)) that is valid in the moment-space interval 

1 


1 < X = \nn < t = 


This perturbative asymptotic representation is 


2aE 


( 10 ) 


N{t)+1 p+1 


Eij{x, a) ~ Eij{x, a, N{t)) = 2Cij XI 


,x-' 


( 11 ) 


p=i 


j=0 


Here 




( 12 ) 






and r(l+z) = where F is the Euler gamma function. The number of perturbative 

terms N{t) in Eq. (pd^ is obtained by optimizing the asymptotic approximation 


E{x, a) — E{x, a, N(t)) 


= minimum. 


(13) 


Because of the range of validity in Eq. m and owing to the optimization Eq. (0 , terms 
in the exponent of the form In^ n are of order unity, and terms with fewer powers of 
logarithms, a^\n^~^n, are negligible. The optimization assures us that the coefficients 
of the various terms are benign. Resummation is completed in a hnite number of steps. 
With a two-loop expression for the running coupling strength, all monomials of the form 
n, In^ n are produced in the exponent of Eq. ©• Because of the restricted 
leading-logarithm universality between the tt and ll processes, we discard monomials of the 
form n in the exponent. 

The exponent we use is the truncation 

N{t)+l 

Eij{x,a,N) = 2Cij (14) 

p=i 

with the coefficients Sp = Sp+i^p = b2~^2P/p{p + 1). The number of perturbative terms N{t) 
is a function of only the top quark mass m. This expression contains no factorially-growing 
(renormalon) terms. The perturbative region of phase space is far removed from the part of 
phase space in which renormalons could be influential. 

In Fig. 2 we illustrate the validity of the asymptotic approximation for a value of t 
corresponding to m = 175 GeV. Optimization works perfectly, with N{t) = 6, and the 
plot demonstrates the typical breakdown of the asymptotic approximation if N is allowed 
to increase beyond N{t). This rise represents the exponential growth of the infrared (IR) 
renormalons, the (p — 1)! growth in the second term of Eq. (^). As long as n is in the 
interval of Eq. ([I0|), all the members of the family in n are optimized at the same N{t), 
showing that the optimum number of perturbative terms is a function of f, i.e., of m only. 

It is valuable to stress that we can derive the perturbative expressions, Eqs. (0), (0. 
and (0)> from the unregulated exponent Eq. (|^) without the PVR prescription, although 
with less certitude. We discuss this point in some detail in Sec. III.B of our long paper [§. 
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After inversion of the Mellin transform from moment space to the physically relevant 
momentum space, the resummed partonic cross sections, including all large threshold cor¬ 
rections, can be written in the form of Eq. (|^), but with the hard kernel replaced by the 
resummed form 

H^j{z,a)= / ^ Qj{x, a) . (15) 

"'o i=o 

The leading large threshold corrections are contained in the exponent Eij{x,a), a calcu¬ 
lable polynomial in x. The functions {Qj{x,a)} arise from the analytical inversion of 
the Mellin transform from moment space to momentum space. These functions are pro¬ 
duced by the resummation and are expressed in terms of successive derivatives of E\ 
Pk{x,a) = d^E{x,a)/k\d^x. Each Qj contains j more powers of a than of x so that 
Eq. (pTSD embodies a natural power-counting of threshold logarithms. 

The functional form of Eij for tt production is identical to that for ll production, except 
for the identihcation of the two separate channels, denoted by the subscript ij. However, 
only the leading threshold corrections are universal. Final-state gluon radiation as well 
as initial-state/hnal-state interference effects produce subleading logarithmic contributions 
that differ for processes with different hnal states. Accordingly, there is no physical basis for 
accepting the validity of the particular subleading terms that appear in Eq. (|T^. Among 
all {Qj} in Eq. (|^, only the very leading one is universal. This is the linear term in Pi 
contained in Qo, that turns out to be Pi itself. Since we intend to resum only the universal 
leading logarithms, we retain only Pi. Hence, Eq. (pAl) can be integrated explicitly, and the 
resummed version of Eq. is 


a. 


R;pert 


_ C^max n /I / 1 \ \ ^ 

{ri,m^) = / z). 


( 16 ) 


We have inserted an upper limit of integration, Zmax, in Eq. (K). This upper limit is 
set by the boundary between the perturbative and non-perturbative regimes. An intuitive 
dehnition of the perturbative region, where inverse power terms are unimportant, is provided 
by the inequality < 1- This inequality is identical to the expression in moment space. 
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Eq. ( 0 ). with the identihcation n = In momentum space, the same condition is 

realized by the constraint that all {Qj}, J > 1 be small compared to Qq. From the explicit 
expressions 0 for the {Qj}, one may show that this constraint corresponds to 

Equation ([T7|) is equivalent to the requirement that terms that are subleading according to 
perturbative power-counting are indeed subleading numerically; Eq. (H) is the essence of 
perturbation theory in this context. It assures us that our integration is carried out only 
over a range in which poorly specihed subleading terms would not contribute signihcantly 
even if they had been retained. 

As remarked above, we accept only the perturbative content of principal-value resum¬ 
mation, and our cross section is evaluated accordingly. Specihcally, we use Eq. (^) with 
the upper limit of integration, Zmax, calculated from Eq. (p!7[) . The upshot is an effective 
threshold boundary on the integral over the scaled subenergy variable rj, but one that is cal¬ 
culable, not arbitrary. While reminiscent perhaps of the cutoff used in the LSvN approach, 
our threshold boundary has a very different and well defined origin. Our perturbative re¬ 
summation probes the threshold down to the point rj > rjQ = {1 — Zmax)/‘^- Below this 
value, perturbation theory, resummed or otherwise, is not to be trusted. For a top mass m 
= 175 GeV, we determine that the perturbative regime is restricted to rj > 0.007 for the 
qq channel and rj > 0.05 for the gg channel. These numbers may be converted to more 
readily understood values of the subenergy above which we judge our perturbative approach 
is valid: at m = 175 GeV, these are 1.22 GeV above the threshold in the qq channel, and 
8.64 GeV above the threshold in the gg channel. The difference reflects the larger color 
factor in the gg case. A larger color factor makes the non-perturbative region larger. (One 
could attempt to apply Eq. ([T6|) all the way to z^ax = 1, he., to g = 0, but one would then 
be using a model for non-perturbative effects, the one suggested by PVR, below the region 
justified by perturbation theory.) We note that the value 1.22 GeV in the qq channel is 
comparable to the decay width of the top quark, T(t —^ bW~^) = 1.55 GeV. The width itself 
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provides a natural definition of the minimum non-perturbative region. The two independent 


determinations of the non-perturbative region are in agreement |11 


IV. PHYSICAL CROSS SECTION 


In order to achieve the best accuracy available we wish to include in our predictions as 
much as is known theoretically. Our hnal resummed partonic cross section can therefore be 
written 


/i^) = /i^) - {r], 


+ ■ ( 18 ) 


R\pert 


The second term is the part of the partonic cross section up to one-loop that is included 
in the resummation, while the last term is the exact one-loop cross section [Q. To obtain 
physical cross sections, we insert Eq. (p!8D into Eq. (0), and we integrate over rj. Other 
than the heavy quark mass, the only undetermined scales are the QCD factorization and 
renormalization scales. We adopt a common value fi for both, and we vary this scale over 
the interval /i/m G {0.5, 2} in order to evaluate the theoretical uncertainty of the numerical 


predictions. We use the CTEQ3M parton densities [|^ 


A quantity of phenomenological interest is the differential cross section 
integral over r] is the total cross section. In Fig. 3 we plot these distributions for m = 175 
GeV, y/S = 1.8 TeV, and = m. The full range of rj extends to 25, but we display the 
behavior only in the near-threshold region where resummation is important. We observe 
that, at the energy of the Tevatron, resummation is signihcant for the qq channel and less so 
for the gg channel. In Fig. 1, the dotted curve shows that our final resummed cross section 
in the qq channel, after integration over all r/, lies about half-way between the cross sections 
obtained from the near-threshold leading logarithms at orders 0{a^) and The latter 

have been integrated over the region Q < z <1. 

We display our inclusive total production cross section as a function of the heavy quark 
mass in Fig. 4. The central value of our predictions is defined as the value obtained with 
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the choice /i/m = 1, and the lower and npper limits are the maximnm and minimnm of the 
cross section in the range of the hard scale /i/m G {0.5,2}. This dehnition of the central 
valne is common, bnt it results here in an asymmetric uncertainty estimate; the extent of the 
range above the central value is smaller than that below. At m = 175 GeV, the full width 
of the uncertainty band is about 10% . In Fig. 5, we show the variation of our resummed 
cross section as the value of the renormalization/factorization scale /i is changed. As is 
to be expected, less variation with /i is evident in the resummed cross section than in the 
next-to-leading order cross section, also shown in Fig. 5. We remark that the cross section 
reaches its maximum at a value of fi just slightly larger than m/2. We consider that the 
variation of the cross section over the range /r/m G (0.5, 2} provides a good overall estimate 
of uncertainty. For comparison, we note that over the same range of /i, the strong coupling 
strength a varies by ±10% at m = 175 GeV. Using a different choice of parton densities 
, we hnd a 4% difference in the central value of our prediction p[ at m = 175 GeV. 
A comparison of the predictions |Q| in the MS and DIS factorization schemes also shows a 
modest difference at the level of 4%. 

In estimating uncertainties, we do not consider explicit variations of our non-perturbative 
boundary, expressed through Eq. (0- For a fixed m and /x, Eq. (0) is obtained by enforcing 
dominance of the leading hard kernel (determined through perturbative power-counting) over 
the subleading ones, all of which are calculable. Therefore, Eq. ([T7|) is derived and is not a 
source of uncertainty. However, at hxed m, the boundary necessarily varies as /i and thus a 
vary. 

Our prediction of the cross section in Fig. 4 is in agreement with the data on top quark 
production [^j. We hnd (j^*(m = 175 GeV, VS = 1.8 TeV) = 5.52lo;42 pb. The central value 
of this cross section is larger than the next-to-leading order value at /x = m by about 9%. 

Extending our calculation at y/S = 1.8 TeV to much larger values of m than shown in 
Fig. 4, we hnd that resummation in the principal qq channel produces enhancements over the 
next-to-leading order cross section of 21%, 26%, and 34%, respectively, at m = 500, 600, and 
700 GeV. The reason for the increase of the enhancements with mass at hxed energy is that 
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the threshold region becomes increasingly dominant. Since the qq channel also dominates 
in the production of hadronic jets at very large values of transverse momenta, we suggest 
that on the order of 20% of the excess cross section reported by the CDF collaboration 
may well be accounted for by resummation. 

The top quark cross section increases quickly with the energy of the pp collider. We 
provide predictions in Fig. 6 for an upgraded Tevatron operating at y/S = 2 TeV. We 
determine a^^{m = 175 GeV, = 2 TeV) = 7.56lo;55 pb. The 2 pb increase in the 
predicted top quark cross section over its value at y/S =1.8 TeV is about a 37% gain. The 
central value rises to 22.4 pb at y/S = 3 TeV and 46 pb at y/S = 4 TeV. For a hxed mass 
of the heavy quark, the fraction of the cross section supplied by the gg subprocess increases 
rapidly. For m = 175 GeV, this fraction is about 15% at y/S = 2 TeV and 51% at y/S = 4 
TeV. 

Turning to pp scattering at the energies of the Large Hadron Gollider (LHG) at GERN, 
we note a few signihcant differences from pp scattering at the energy of the Tevatron. The 
dominance of the qq production channel is replaced by gg dominance at the LHG. Owing to 
the much larger value of y/S, the near-threshold region in the subenergy variable is relatively 
less important, reducing the signihcance of initial-state soft gluon radiation. Lastly, physics 
in the region of large subenergy, where straightforward next-to-leading order QGD is also 
inadequate, becomes signihcant for tt production at LHG energies. Using the approach 
described in this paper, we estimate a^^{m = 175 GeV, y/S = 14 TeV) = 760 pb. 


V. OTHER METHODS OF RESUMMATION 


The groups of Laenen, Smith, and van Neerven (LSvN) and of Gatani, Mangano, Nason, 
and Trentadue (GMNT) have also published predictions for the total cross section based on 
resummation of initial state soft gluon radiation. At m = 175 GeV and y/S = 1.8 TeV, 
the three values are: (j**(BG 0J^) = 5.52to;42 pb; (j**(LSvN |]^) = 4.95lo;4o pb; and 
cr“(GMNT §) = 4.75l)];i pb. From the purely numerical point of view, all three pre- 


14 





dictions agree within their estimates of theoretical uncertainty. However, the resummation 
methods differ, the methods for estimating the uncertainties differ, and different parton sets 
are used. Comparing with LSvN |^, we hnd that our central values are 10 —14% larger, and 
our estimated theoretical uncertainty is 9 — 10% compared with their 28% —20%. The larger 
central value is attributable, in part, to the use of different parton densities; our Born cross 
section is about 3 — 5% larger than the LSvN Born cross section. However, it is the choice of 
the infrared cutoff po in the LSvN method that controls the size of their cross section. The 
cutoff Ho is selected so that the resummed cross section is about equal to the next-to-next- 
to-leading order leading-logarithm cross section obtained from Eq. (^). In contrast, 

in our approach the non-perturbative boundary z^ax is derived within the context of the 
calculation by the requirement that the universal leading-logarithmic terms be dominant. 
There is no a priori reason that our resummed result should be only 10% greater than the 
next-to-leading order cross section at m = 175 GeV and y/S = 1.8 TeV. As such, we regard 
the approximate agreement of our result and that of LSvN as somewhat fortuitous. Both 
the central value and the band of uncertainty of the LSvN predictions are sensitive to their 
infrared cutoffs, as we described previously |0]. 

From a theoretical point of view, study of the variation of the predicted cross section with 
the hard scale /x, illustrated here in Fig. 5, is important because it reflects the stability of the 
calculation under changes of a perturbative but not directly determinable renormalization- 
factorization scale. One of the advantages of a resummation calculation should be diminished 
dependence of the cross section on /x, less variation than is present in a hxed-order calcu¬ 
lation. To estimate theoretical uncertainty, we use the standard /x variation, and we hnd a 
band of uncertainty of about 10% at m = 175 GeV and y/S = 1.8 TeV. The LSvN group 
obtain their uncertainty primarily from variations of their infrared cutoff whose role is to 
measure ignorance of non-perturbative effects in that approach. 

The group of Gatani, Mangano, Nason, and Trentadue (GMNT) calculate a central 
value of the resummed cross section (also with /x/m = 1) that is less than 1% above the exact 
next-to-leading order value. There are similarities and differences between our approach to 
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resummation and the method of Ref. [^]. We both begin in moment space with the same 
universal leading-logarithm expression, but differences occur after the transformation to 
momentum space. In this paper, we set aside comments on mathematical aspects of their 
procedure and focus instead on phenomenological issues of interest. As remarked above, the 
Mellin transformation generates subleading terms in momentum space. The suppression of 
the effects of resummation arises from the retention in Ref. of numerically signihcant 
non-universal subleading logarithmic terms. 

CMNT choose to retain all of these inasmuch as they perform the Mellin inversion numer¬ 
ically. Instead, in keeping with the fact that subleading logarithmic terms are not universal, 
we retain only the universal leading logarithm terms in momentum space, and we restrict 
our phase space integration to the region in which the subleading terms would not be nu¬ 
merically signihcant regardless. The differences in the two approaches can be stated more 
explicitly if we examine the perturbative expansion of the kernel T-C^{z,a), Eq. (p!5[) . If, 
instead of restricting the resummation to the universal leading logarithms only, we were to 
use the full content of Eq. (0, we would arrive at an analytic expression that is equivalent 
to the numerical inversion of Ref. ^ , 


Tfg ~ 1 + 2aC, 




ln^(l - z) + 27 £;ln(l - z) 


+ 0{a^ 


(19) 


where 'Je is Euler’s constant. In terms of this expansion, in our work we retain only the 
leading term ln^(l — z) at order a, but CMNT retain both this term and the subleading 
term 27 £;ln(l — z). Indeed, if the subleading term 27 £;ln(l — z) is discarded in Eq. (|T^, 
the residuals dehned in Ref. increase from 0.18% to 1.3% in the qq production 

channel and from 5.4% to 20.2% in the gg channel [^. After addition of the two channels, 
the total residual grows from the negligible value of about 0.8% cited in Ref. 

to the value 3.5%. While still smaller than the increase of about 9% that we obtain, the 
increase of 3.5% vs. 0.8% shows the substantial influence of the subleading logarithmic terms 
retained in Ref. |^. 

We judge that it is preferable to integrate over only the region of phase space in which the 
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subleading term is suppressed numerically. Our reasons include the fact that the subleading 
term is not universal, is not the same as the subleading term in the exact 0{a^) calculation, 
and can be changed if one elects to keep non-leading terms in moment space. The subleading 
term is negative and numerically very signihcant when it is integrated throughout phase 
space (i.e., into the region of z above our Zmax)- In the qq channel at m = 175 GeV and 
y/S = 1.8 TeV, its inclusion eliminates more than half of the contribution from the leading 
term. In our view, the presence of numerically signihcant subleading contributions begs 
the question of consistency. A further justihcation for the retention of only the leading 
term is that it approximates the exact next-to-leading order result well, as shown in Fig. 1. 
The choice made in Ref. [Q reproduces only one-third of the exact next-to-leading order 
result. The inhuence of sub leading terms is amplihed at higher orders where additional 
subleading structures occur in the approach of Ref. 0] with signihcant numerical coefficients 
proportional to vr^, C(3)) cind so forth. We present a more detailed discussion of this issue 
in the next section. 

VI. FURTHER DISCUSSION OF THE CMNT APPROACH 

In this section we oher a more systematic analysis of the role played in the approach 
of Ref. @] by non-universal subleading logarithms. We are interested in expansions of the 
resummed momentum-space kernel, Eq. ([I5|) , up to two loops. Therefore, the corresponding 
cross sections are integrable down to threshold, Zmax = 1 and rj = 0. As we will see, though, 
the ehects of the various classes of logarithms are pronounced if one continues the region of 
integration outside our perturbative regime. 

In moment space, the exponent to two-loops is obtained from Eq. (O): 

a) = ga{s 2 ^ix‘^ + Siqx - 1 - So,i) + ga‘^{s3^2X^ + ^ 2 ^ 2 ^'^ + 51 , 23 ^ + 50,2), ( 20 ) 

with g = 2Cij and x = Inn. The corresponding hard kernels in momentum space can 
be derived from Eq. (|^), according to the formulas (91) through (94) of Ref. [Q]. In the 
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notation of Ref. [Q], we retain terms up to those that are linear in P 2 . Alternatively, one can 
perform the analytical Mellin inversion directly, beginning with Eq. (0)' The two methods 
provide identical results down to the monomials xla^; Xz = ln(l/(l — z)). Here we quote 
results based on the explicit inversion of Eq. (pOD. (The monomial can be obtained 
also in the hrst approach if we keep the quadratic term in P 2 .) After a trivial integration is 
performed, the results for the one- and two-loop hard kernels are 

= xla{gs2,i} + Xza{g{si^i + 2ciS2,i)} , (21) 

and 


Tfd) = xla^{g‘^sl^^/2} + + g‘^{s2,iSi^i + 2ci4i)} 

+x\a^{g{s2^2 + 3ciS3^2) + 5'^('5iq/2 -l- 3ciSi^iS2,i -l- S2,iSo,i -l- S2^^[6c2 — vr^]} 
+Xz 0 i‘^{g{s \^2 + 2CiS 2,2 + Ss,2(602 — vr^]) 

+fi'^(so,lSi_l -|- 2 ciSo,iS2,1 + + S2,iSi^i[ 6C2 — TT^] -f S2,l[l2C3 — 27r^Ci])} . (22) 


All the constants are dehned in Eqs. (^) and (|T^. We remark that Eq. 0) includes a 
leading logarithmic term, x^a, as well as a next-to-leading term, x^a. 

The question we now address is whether it is justihed and meaningful to retain all of the 
terms in Eqs. ( [2T|) and (P^ in the computation of the resummed cross section. The issue 
has to do with what one intends by resummation of leading logarithms. We use the term 
leading logarithm resummation to denote the case in which the moment space exponent, 
Eq. (p0|), contains only the constants E^l = {'Sp+i,p,0}. This is also what is done in the 
method of Ref. @], and the exponent in moment space in their work is identical to that used 
for our predictions, Eq. ( 0 )- However, in contrast to our expression in momentum space. 


Eq. (IT), the corresponding expression in momentum space of Ref. [Q includes the numerical 
equivalent of all terms in Eqs. (0 and (1^) that are proportional to Sp+i^p. 

If expressed analytically, the corresponding “LL” hard kernels in the method of Ref. 0 


are 


= ^1(^9 - Xza2g'jE, 


(23) 
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and 


= xla^g^l2 + xla^{2gb2/?> - 2^Eg‘^} 

+xla^{-2gb2iE + /[StI - 7rV2]} 

■^x^a^{2gb2[^'yE “ 7rV2]/3 + - 27 ^ - 4C(3)]} , (24) 

where ({s) is the Riemann zeta function; (^(3) = 1.2020569. Evaluating the expressions 
numerically for the qq channel, we obtain 

= xla X 2.66666 - x^a x 3.07848, (25) 


and 


fdll = xW X 3.56655 - xW x 4.80189 - xW x 33.88456 - x.o? x 9.82479 


(26) 


Apart from the leading monomials that are the same as those in our approach, Eqs. 
and ( p 6 |) include a series of subleading terms, each of which has a signihcant negative coef- 
hcient. In practice, these subleading terms in the approach of Ref. suppress the effects 
of resummation essentially completely. One of the effects of this suppression is that the 
resummed partonic cross section is smaller ihsu its next-to-leading order counterpart in the 
neighborhood of 7 = 0 . 1 , a region in which the next-to-leading order partonic cross section 
takes on its largest values. This point is illustrated in Fig. 3 of the second paper in Ref. [Q. 

Although the specihc set of subleading terms in Eqs. (^) and (|^) is generated in the 
inversion of the Mellin transform, we would argue that the terms are accidental, at best. 
Our reasoning is based on an examination of the exact next-to-leading order calculation 
of the cross section for heavy quark production and of similar calculations of the Drell-Yan 
process up to two-loops. First, terms involving '^e do not appear in the exact next-to-leading 
order calculation of the hard part, since they are removed in the specihcation of the MS 
factorization scheme. Therefore, the term proportional to 7 ^ in Eq. (^) is suspect. Second, 
if we extract the specihc value of the subleading logarithm from the full 0{a^) next-to- 
leading order calculation P], we hnd |T^ Xza{2g — 4:l/6) instead of the term —Xz(y2g'yE in the 
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equivalent CMNT Eq. ([^). Referring to Eq. (pSQ , we remark that instead of the numerical 
coefficient 3.07848, one would have the smaller value 1.5 if the subleading logarithm of the 
exact 0{a^) calculation were used. Thus, not only is the 0(a) subleading term retained in 
the approach of Ref. ^ different from that of the exact calculation, it is numerically about 
twice as large. Third, we would claim that the results of a LL resummation should not rely 
on the subleading structures in any signihcant way. However, in the approach of Ref. 0, 
Eq. which is the one-loop projection of their resummed prediction, reproduces only 1/3 

of the exact 0{a^) enhancement, the other 2/3 being cancelled by the second (non-universal) 
term of Eq. (^3]). Correspondingly, the method of Ref. 0 fails an important consistency 
check: it sets out to resum the threshold corrections responsible for the large enhancement 
of the cross section at next-to-leading order; in the end, it does not reproduce most of this 
enhancement. 


Addressing questions associated with the y® terms CMNT examine a type of NLL 
resummation in the second paper of Ref. [Q. In this NNL resummation, the {sp+i^p, Sp^p} 
terms are retained in the exponent of Eq. (P0[). The corresponding hard kernels become 

(27) 


WjJh = x^ag, 


and 

'M'nll = xta'^g ''‘/2 + 1^^2962/3 - + irV2] - + >r^/ 3 ] + gnc( 3 )} . 

(28) 


Comparing Eqs. (pTI) and with Eq. (^, we observe that Eq. (|27|) is identical to the one- 
loop projection of our hard kernel. As shown in Ref. it yields an excellent approximation 
to the exact next-to-leading order cross section. On the other hand, our two-loop projection 
contains only the hrst two terms of Eq. (^). The term proportional to is present in our 
case, along with the leading term proportional to because it comes from the leading 

logarithms in the exponent E{n), through one-loop running of the coupling strength. In 
contrast to Eq. p^), Eq. (^) relegates the influence of the ambiguous constant coefficients 
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to lower powers of Xz (but with larger negative coefficients). In the amended scheme of 
Ref. the unphysical 7 ^; terms are still present in the two-loop result, Eq. (pSf), along with 
TT^ and C (3) terms that may be expected but whose coefficients have no well dehned physical 
origin. Recast in numerical form, Eqs. (|27|) and (p 8 D become 


'^ n\l = X 2.66666, 


(29) 


and 


wSh = X 3.55555 + x 3.40739 - xia^ x 37.46119 - x 54.41253 


(30) 


We call attention to the significant difference between the coefficients of all but the very 
leading power of Xz in Eqs. (p5[) and (|26|) with respect to those in Eqs. (|^) and (^), and 


to the fact that the numerical coefficients grow in magnitude as the power of x^ decreases. 

Using their NLL amendment, CMNT End that the central value of their resummed cross 
section exceeds the next-to-leading order result by 3.5% (both qq and gg channels added). 
This increase is about 4 times larger than the central value of the increase obtained in their 
first method, closer to our increase of about 9%. The reason for the significant change of 
the increase resides with the subleading structures, viz., in the differences between the LL 
version Eqs. (p5|) and (|^) and the NLL version Eqs. (|^) and (^). The subleading terms at 
two-loops cause a total suppression of the two-loop contribution (in fact, that contribution is 
negative), if one integrates all the way into what we call the non-perturbative regime. This 
suppression explains why an enhancement of only 3.5% is obtained in the amended method 
of Ref. 0], rather than our 9%. 

CMNT argue that retention of their subleading terms in momentum space is important 
for “energy conservation”. By this statement, they mean that one begins the formulation of 
resummation with an expression in momentum space containing a delta function represent¬ 
ing conservation of the fractional partonic momenta. In moment space, this delta function 
subsequently unconvolves the resummation. Therefore, when one inverts the Mellin trans¬ 
form to return to momentum space, the full set of logarithms generated by this inversion are 
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required by the original energy conservation. This line of reasoning would be compelling if 
the complete exponent E{n) in moment space were known exactly, i.e., if the resummation 
in moment space were exact in representing the cross section to all orders. However, the 
exponent is truncated in all approaches, and knowledge of the logarithms it resums reliably 
is limited both in moment and in momentum space. Hence, the set of logarithms produced 
by the Mellin inversion in momentum space should also be restricted. In our approach en¬ 
ergy conservation is obeyed in momentum space consistently with the class of logarithms 
resummed. On the other hand, in the method of Ref. @] , knowledge is claimed of all loga¬ 
rithms generated from the Mellin inversion, despite the fact that the truncation in moment 
space makes energy conservation a constraint restricted to the class of logarithms that is 
resummable, i.e., a constraint restricted by the truncation of the exponent E{n). The two 
approaches would be equivalent provided a constraint be in place on the effects of subleading 
logarithms. This constraint is precisely our Eq. (0)- By contrast, no such constraint is fur¬ 
nished in Ref. IQ . For this reason the results of Ref. ^ are numerically unstable if one set of 
the logarithms generated in momentum space is adopted as “the set corresponding to energy 
conservation”, and then compared with another set, produced by a different truncation of 
E{n). 

We have identified the terms responsible for the difference between our answer for the 
resummed cross section and that of Ref. [Q. These differences reside with subleading log¬ 
arithms whose presence is not substantiated by physical arguments. The essence of our 
determination of the perturbative regime in Eq. (pTl) is precisely that, in this regime, sub¬ 
leading structures are also numerically subleading, whether or not the classes of subleading 
logarithms coming from different truncation of the master formula for the resummed hard 
kernel, Eq. (p!5D, are included. The results presented in Fig. 11 of Ref. 0, show that if we 
alter our resummed hard kernel to account for subleading structures but still stay within our 
perturbative regime, the resulting cross section is reduced by about 4%, within our band of 
perturbative uncertainty. 

A criticism is that of putative “spurious factorial growth” of our resummed cross 
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section, above and beyond the infrared renormalons that are eliminated from our approach. 
The issue, as we demonstrated in Eq. (29) of Ref. 0, can be addressed most easily if we 
substitute any monomial appearing in Eq. (p2|), symbolically a'^c{l,m)\vL Xz, into Eq. 
and integrate over 2 ;: 

a"^c{l,m) f dzW Xz = - Zmin)) ■ (31) 

^min j — Q 

For the purposes of this demonstration we set dh = 1. The coefficients c{l,m) can be read 
directly from Eq. (^^ . For the leading logarithmic terms, 

c(2m, m) oc 1/m!, (32) 


where this factorial comes directly from exponentiation. After the integration over the entire 
z-range, the power of the logarithm in Xz becomes a factorial multiplicative factor, /!. The 
presence of l\ follows directly from the existence of the powers of Ina;^ that are present 
explicitly in the finite-order result in pQCD and is therefore inevitable. If this exercise 
is repeated, but with the range of integration in Eq. (|31|) constrained to our perturbative 
regime, one obtains the difference between the right-hand-side of Eq. m and a similar 
expression containing Zmax- The result is numerically smaller, but both of the pieces are 
multiplied by /!. 

The factorial coefficient l\ is neither the only nor the most important source of en¬ 
hancement. For the leading logarithms at two-loop order, I = 2m = 4, and the overall 
combinatorial coefficient from Eqs. (|3T|) and (|3^) is (2m)!/ml = 12. For comparison, at 
representative values of rj near threshold, rj = 0.1 and 0.01, the sum of logarithmic terms 
in Eq. (^T]) provides factors of 16.1 and 314.3, respectively. Similarly, the (multiplicative) 
color factors at this order of perturbation theory are {2CijY = 7.1 and 36 for the qq and 
gg channels, respectively. All of these features are connected to the way threshold loga¬ 
rithmic contributions appear in hnite-order pQCD and how they signal the presence of the 
non-perturbative regime. Thus, preoccupation with the l\ factor seems misplaced. 

The phrase “spurious factorial growth” appears to rename the logarithmic enhancements 
present in Eqs. (^p and (|^) , after the integral over 2 ;. On the other hand, according to our 
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understanding, the claim |Q of “absence of factorial growth” is based on the use in Ref. 0] 
of Eq. (p6D for their main predictions, an expression that contains non-universal sublead¬ 
ing logarithms, all with signihcant negative coefficients. Mathematically, factorial growth 
is present for each of the powers of the logarithm in Eq. (|^, since these monomials are 
irreducible (linearly independent). Absence of factorial growth based on a numerical cancel¬ 
lation between various classes of logarithms, most of them with physically unsubstantiated 
coefficients, appears to us to be an incorrect use of terminology, rather than a transparent 
expression of the mathematics. 

From a purely phenomenological point of, one cannot claim that a 9% increase of the top 
quark cross section at m = 175 GeV and y/S = 1.8 TeV reveals factorial growth but that 
an 0.8% increase does not. In the approach taken in Ref. [Q, the effects of resnmmation 
are suppressed by a series of subleading logarithms with large negative coefficients. If there 
is no physical basis for preference of Eqs. (|^) and (|^) over Eqs. (^71) and (|^) , as the 
authors of Ref. [Q seem to suggest, then the difference in the resulting cross sections can 
be interpreted as a measure of theoretical uncertainty. This interpretation would not justify 
hrm conclusions of a minimal 0.8% increment based on the choice of Eqs. (p3D and (p^). 

As remarked in Sec. IV, the value quoted in Ref. for the physical cross section at m = 
175 GeV and \/S = 1.8 TeV, including theoretical uncertainty, lies within our uncertainty 
band. Therefore, the numerical differences between us for the specihc case of top quark 
production at the Tevatron have little practical signihcance. However, there are important 
differences of principal in our treatment of subleading contribntions that will have more 
signihcant conseqnences for predictions in other processes or at other valnes of top mass 
and/or at other energies, particularly in reactions dominated by gg subprocesses. 


VII. DISCUSSION AND CONCLUSIONS 


In this paper, we present and discuss the calculation of the inclusive cross section for 
top qnark production in perturbative QGD, including the resummation of initial-state gluon 
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radiation to all orders in the strong coupling strength. The advantages of the perturba¬ 
tive resummation method we espouse are that there are no arbitrary infrared cutoffs 
and there is a well-dehned perturbative region of applicability where subleading logarithmic 
terms are numerically suppressed. Our theoretical analysis shows that perturbative resum¬ 
mation without a model for non-perturbative behavior is both possible and advantageous. 
In perturbative resummation, the perturbative region of phase space is separated cleanly 
from the region of non-perturbative behavior. 

When evaluated for top quark production at y/S =1.8 TeV, our resummed cross sections 
are about 9% above the next-to-leading order cross sections computed with the same parton 
distributions. The renormalization/factorization scale dependence of our cross section is 
fairly flat, resulting in a 9 — 10% theoretical uncertainty. This variation is smaller than 
the corresponding dependence of the next-to-leading cross section, as should be expected. 
Our perturbative boundary of 1.22 GeV above the threshold in the dominant qq channel is 
comparable to the hadronic width of the top quark, a natural dehnition of the perturbative 
boundary. Neither this, nor the somewhat larger value of 8.64 GeV above threshold in the 
gg channel, associated with the larger color factor in the gg channel, is “unphysically large” 
In recent papers |^, the authors state that the increase in cross section they hnd with 
their resummation method is no more than 1% over next-to-leading order. The numerical 
difference in the two approaches boils down to the treatment of the subleading logarithms, 
which can easily shift the results by a few percent, if proper care is not taken. Our approach 
includes the universal leading logarithms only while theirs includes non-universal subleading 
structures which produce the suppression they hnd. In Sec. VI, we explain why we judge 
that that our treatment of the subleading structures is preferable. 

Our theoretical analysis and the stability of our cross sections under qi variation provide 
conhdence that our perturbative resummation procedure yields an accurate calculation of the 
inclusive top quark cross section at Tevatron energies and exhausts present understanding of 
the perturbative content of the theory. Our resummed top quark cross section is about 9% 
above the next-to-leading order cross section with an estimated theoretical uncertainty of 


25 


9 — 10%, associated with /i variation. An entirely different procedure to estimate the overall 
theoretical uncertainty is to compare our enhancement of the cross section above the next- 
to-leading order value to that of Ref. [Q, again yielding about 10%. An interesting question 
is whether theory can aspire to an accuracy of better than 10% for the calculation of the top 
quark cross section. To this end, a more complete mastery of subleading logarithms would 
be desirable, perhaps requiring a formidable complete calculation at next-to-next-to-leading 
order of heavy quark production, to establish the possible pattern of subleading logarithms, 
and resummation of both leading and subleading logarithms . 

Our prediction agrees with data, within the large experimental uncertainties. Despite 
the different treatment of subleading terms, our calculation of the inclusive cross section 
for top quark production at the Fermilab Tevatron and that of Ref. ^ fall within the 
estimated uncertainties of each other. If a cross section signihcantly different from ours is 
measured in future experiments at the Tevatron with greater statistical precision, we would 
look for explanations in effects beyond QCD perturbation theory. These explanations might 
include unexpectedly substantial non-perturbative effects or new production mechanisms. 
An examination of the distribution in rj might be revealing. 

In this paper, we concentrate on the all-orders summation of large logarithmic terms that 
are important in the near-threshold region of small values of the scaled partonic subenergy, 
?7 —> 0. Our specihc case is top quark production at the Fermilab Tevatron collider. Other 
processes for which threshold resummation and our methodology will also be pertinent in¬ 
clude the production of hadronic jets that carry large values of transverse momentum and 
the production of pairs of supersymmetric particles with large mass. There is a comple¬ 
mentary region of large 77,77 —>■ 1, in which the resummation of different large logarithms 
may also be important. The production of heavy quarks Q in the limit that the hadronic 
center-of-mass energy is much larger than the quark mass provides an example. The domi¬ 
nant production channel is gg —> QQX ; the ratio of the next-to-leading order partonic cross 
section divided by its leading-order approximation is very large at large 77. Correspondingly, 
the fixed-order cross section will not offer a reliable prediction, and an all-orders approach is 
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called for Particular cases include the total cross sections for bottom quark production 
at the Tevatron and top quark production at the CERN Large Hadron Collider. 
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FIGURES 

FIG. 1. Physical cross sections in the qq channel as a function of the heavy quark mass, in the 
MS scheme. The solid lines denote the hnite-order partial sums of the universal leading-logarithmic 
contributions from the explicit 0{a^) and 0{a‘^) calculations for the tt and Drell-Yan processes, 
respectively. Lower solid: middle solid; upper solid: The dashed curve 

represents the exact next-to-leading order calculation for ti production, in excellent agreement 
with The dotted curve is our resummed prediction. 

FIG. 2. Optimum number of perturbative terms in the exponent with PVR. The solid family 
is for PVR and the dashed set for the perturbative approximation, both families increasing, for 
parametric values n = 10, 20,30,40. 

FIG. 3. Differential cross sections da/dri for pp tiX at \/5 = 1.8 TeV and m =175 GeV 
in the MS-scheme for (a) the qq and (b) the gg channel; Born (dotted), next-to-leading order 
(dashed), and resummed (solid). 

FIG. 4. Inclusive cross section for heavy quark production at ^/S =1.8 TeV in the MS scheme. 
The dashed curves show our perturbative uncertainty band, while the solid curve is our central 
prediction. 

FIG. 5. Renormalization/factorization hard scale dependence of the resummed (solid) and 
next-to-leading order (dashed) cross sections at ^/S = 1.8 TeV for m =175 GeV. 

FIG. 6. Inclusive cross section for heavy quark production at ^/S = 2 TeV in the MS scheme. 
The dashed curves show our perturbative uncertainty band, while the solid curve is our central 
prediction. 
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